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Abstract 

^ i We study the renormalization of the 0(N) model using the auxiliary field formalism (Hubbard- 

' Stratonovich transformation) in the 1/N expansion at finite temperature. We provide the gen- 

■ eral strategy of renormalization for arbitrary order, and make calculation up to next-to- leading 
order. We show that renormalization is possible for any values of the condensates, prove the 
temperature independence of the counterterms and determine the cutoff dependence of the first 
nontrivial counterterm parts. 

D 

r^, - 1 Introduction 

^ I Renormalization of approximations in quantum field theories which involve resummation represents 
Q ■ a complicated, much studied subject. Depending on the system, various techniques are developed. 
O . Renormalization of mass and coupling constant resummation was studied in [1],[2], in the frame- 
work of 2PPI-resummation in |3]. Renormalization of Hartree-Fock resummation is performed in 
[3! , [H] , [H] , [Z] • The task of the renormalization of 2PI resummation was solved by [8j,[9j,[10j, [llj.[12j. 
O , its generalization to n-point irreducible resummations is discussed in [13]. Application of this method 
^ ■ to scalar models [Hj , [IH] , [1^ and gauge theories [I7],[18] was done. The problem of how the 2PI 
^ counterterms can be constructed is discussed in [19]. In [20] a renormalization scheme method was 
proposed to accomplish 2PI resummation and renormalization in a single step. Renormalization of 
rS , Schwinger-Dyson equations and its relation to the 2PI approach, was discussed in [21j.^22j. 

■ Renormalization of 0(N) model beyond the by-now textbook case leading order [23] in 1/A^ 
expansion was studied recently in [21], [2S], using 2PI techniques in [2B]- In [2S] the authors presented 
a technique, with help of which they could provide a renormalized, background dependent free energy 
up to next-to-leading order. Although they used auxiliary field formalism (Hubbard-Stratonovich 
transformation), after a point they switched to a simpler background, eliminating the auxiliary field 
expectation value using its equation of motion (EoM). It is not clear, whether there exists a consistent 
free energy also for the auxiliary field background. The problem is that the auxiliary field represents 
a composite operator of the original fields, and so it may show unusual properties [25].[27j. 
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The primary goal of this paper is to propose a method to consistently determine the auxiliary and 
scalar field dependent renormalized free energy. It is important to define in a physically sensible way 
the pressure (minus the saddle point value of the free energy density) and also the different n-point 
functions at zero external momentum (from the derivatives of the free energy near the saddle point). 
For a consistent renormalization one has to know the (temperature and condensate- independent) 
values of the counterterms in a given regularization. Here we use momentum cutoff regularization, 
and we shall give the leading order values of the counterterms. 

The method proposed here represents an application of the general strategy discussed in [20] . 
We define a generalized renormalized Lagrangian, where also the counterterms can be momentum 
dependent, but the bare theory remains intact. The consistency of such renormalization schemes 
was studied in [20] , here we have to check if the asymptotic behavior of the propagator fulfills the 
generic requirements as derived in [20] . 

The paper is organized as follows. First the definition of the model is presented, then the renormal- 
ized Lagrangian is given in the most general scheme which respects all symmetries of the Lagrangian. 
Then nontrivial background fields are introduced. For renormalization we first fix the renormaliza- 
tion conditions using the classical Lagrangian, then we determine at leading and next-to leading 
order the values of the counterterms. We prove that the counterterms are temperature independent, 
and give formulae for the renormalized free energy. The paper is closed with conclusions. 



2 The renormalized Lagrangian 

The bare Lagrangian of the system reads: 

For tracing the powers, it is advantageous to split the quartic interaction via the Hubbard- 
Stratonovich transformation. The above model is equivalent to 

C ^ lid^md'^i) - ^^1^1 - \f - (2) 
or, in the imaginary time 

where 

In order to be able to do perturbative computations, we should separate in the Lagrangian the 
renormalized and the counterterm parts. We start with the wave function renormalization: 

$ = ^ = Z'JW (5) 

By introducing 

ml = Zrv?, go = ZZ^^'^g, (6) 
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we can write 

Ce = ^{d,md,^d + + ^Xl + ^^0$.$.. (7) 

As next we should renormalize the operators appearing in the Lagrangian. The basic principle 
is that all operators which are scalar according to the symmetry group of the Lagrangian and with 
positive mass dimensions must appear in the renormalized Lagrangian. This means, that we have to 
include, in addition to the operator set of (171), the operators xo and - this latter is necessary, 

since in the actual form of the Lagrangian ([7]) the and the Xo operators are independent, and 

so they must be renormalized independently. It may turn out, that we do not need counterterm for 
all these operators, then we can omit them, but, for sake of completeness, we include all possibility 
here. To allow later a full optimization of the propagator, we introduce a generic kernel and quadratic 
counterterm for the xo fields. With the notations 

Z =1 + 6Z, ml = m^ + 5m^, Z^ = l + 5Z^, (8) 

the renormalized Lagrangian reads 

Ce = i-d^ + m^) + ^xo H{id) Xo + ^^^^^^^ + 

+^iqXo + {-5Zd^ + W) + ^xo 5H{id) Xo + ^Xo$.$. + ^("^^^^O', (9) 

where the first row is the renormalized, the second is the counterterm part. Some iV-power factors 
are introduced just for convenience. The consistency to the previous form of the Lagrangian requires 
(changing to Fourier space for easier writing): 

H{p) + 5H{p) = Z^. (10) 

The choice of the kernel is to large extent arbitrary, since - if the above constraint is fulfilled - any 
choice leaves the bare Lagrangian, and so the physics, untouched. That means that choices which 
depend on the environment (temperature, singlet backgrounds, etc) are also possible. It is important 
to note that, if H{p) is chosen to be 0(N) singlet, the Lagrangian remains 0(N) symmetric, and 
so in any perturbation theory based on this Lagrangian, the consequences of the symmetry (ie. the 
Ward identities) remain valid. 

A different question is whether a momentum-dependent counterterm can result in, at all, a 
consistent perturbation theory. This question was addressed in [20]: if the kernel can be power 
expanded around asymptotic momenta, then the Weinberg theorem is obeyed, so the renormalization 
program can be fully implemented. We will check later that these conditions are really satisfied. 

In this paper we aim to perform expansion with respect to 1/N to first order (ie. up to 0{1/N) 
in the free energy). As it will be later confirmed, we assume the following dependence of the coun- 
terterms on 




OO OO OO 

^'^'^ = T.^^'^l^ ^9 = J2l^^9n, S\=Y.—5Xn (11) 

n=l n=l n=l 
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where 5f is the counterterm for the free energy. Note that in certain cases we have zeroth order 
counterterms as well! 

In the perturbative calculation we fix the splitting of the bare parameters into renormalized 
and counterterm pieces. This separation must be done in a way that the divergences appearing 
in the expressions of physical quantities cancel at each order. This determines the infinite part of 
the counterterms. For the finite parts we require that certain physical quantities have pre-defined 
values (renormalization scheme). For all divergent quantities we should therefore find a physical 
quantity, which fixes the value of the finite part; in our case we have seven potentially divergent 
quantities {q,6f,Sm'^,SZ,SZ^,Sg,SX), so we need seven renormalization conditions. These will be 
defined through the free energy and its derivatives later. 



3 Spontaneous symmetry breaking 

In the ground state of the system the bosonic scalar fields can acquire expectation value. In general 
we can assume 

(Xo) = -iVNX, i^i) = VN^i. (12) 

We will allow only those choices for H{p), which depend on the 0(N) invariant combination of 
the background fields, ie. on J2i ^i^i- Then the complete Lagrangian will be invariant under a 
simultaneous 0(N) rotation of the background and the fluctuation flelds. Then, with an appropriate 
transformation we can achieve that = for i = 1, . . . — 1, while 

i 

As a consequence all results can depend only on J2i ^i^i, in particular the free energy, too. This 
ensures that the flrst and second derivative of the free energy (effective action) for the k,j < N 
modes can be written as f{X, $j) = V{X, $^), and so 

If at the saddle point d<j>V = (we are in the broken phase), the second derivative is a tensor with 
— 1 zero modes. Since the second derivative of the free energy is the inverse propagator at zero 
momentum, we find A^ — 1 zero modes, independently on the value oi X. So the Goldstone theorem 
is satisfied, independently of X. 

We introduce fluctuation flelds, which already have zero expectation value 

Xo = -i\/7VA: + x, <^N^VN<^+g, <^>i = ^i{t = l...N -1). (15) 
We write the resulting Lagrangian as 

r /'class I /'Class, ct , -r lin , /'(S) , pct,! /IcA 

J-E — i-E + '-E + ^E + '-E + '-E ■> UOJ 

where 



£dass _ ^ 



2 2 2 
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/'clasSjCt 

E 



r'Vm 
'-E 

'-E 



$2 \J-X^ + qX + —X^^ H $^ 



mg$ + + + ) ^ + i ( -Xi/(0) + + g - X5i7(0) + ) % 



+ 



5H{id) \x 



+ 



i5g , 2\ 5-^^ / 2\ '^'^ / 2\2 



(17) 



The (^j (z = 1 ... — 1) modes will be called pions. Accordingly, we will denote the tree level 
mass of the pions as: 

ml = m'^ + gX. (18) 

Note that the so-defined pion mass depends on the condensate. The quadratic term C"^ provides 
the propagators, for which the following notations will be used: 
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p"^ + rri^ ig^ 
ig^ H{p) 



in particular 



GxxiP) 



p^ + 



(19) 



(20) 



H{p){p^ + ml)+g'^^^' 

We will compute the constrained free energy or effective potential, which is the IPI effective 
action for constant background; with the above separation we can write 



N -I 



THog + + iTHogdat + < ) + ^ (l - 



IPI 



(21) 



The first term provides the classical result, the rest is due to quantum corrections. The interactions 
are treated through Taylor expansion: 



1 



1 -e" 



IPI 



(22) 



4 Renormalization 



After defining the generic framework we determine the necessary countertcrms. First we fix the 
physics by requiring some renormalization conditions. Then order by order compute the necessary 
diagrams, and require the cancellation of the divergences. 
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4.1 Renormalization conditions 



To fix the renormalization conditions we will require that the radiative corrections do not spoil the 
robust phenomena of the classical theory. To this end we discuss the classical free energy. Without 
resummation it reads 

k = !!^$2 _ 1^2 9^^^2 (23) 

N 2 2 2 ^ ^ 



Its saddle point is at position 

dfci 



dX 



-^min ~l~ 



£ref 2 

In the broken phase, where < and $r 



0, 



dfa 



^ref 

7^ 0, its solution reads 



X„ 



-m 



$2 
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-2m^ 



(24) 



(25) 



The first condition says that the pion mass at the saddle point — + gX^i^ = 0, in accordance 
with the Goldstone theorem. The rho-mass can be determined from the condition that the g-x 



propagator at the saddle point of the free energy has a pole at p 



-m . Its position is found by 



requiring a zero for the determinant of the q-x kernel at the saddle point: 



det 



ig^ H{-ml) 



-miH{-mi)+g^^' 







g'^ 



2 

min 



-2m^ 



(26) 



Here we used that in the classical case H{p) — 1. 

In the symmetric phase, where is positive we obtain X^ 



0, and — m? 



m 



We will require later when we compute radiative corrections, that the position of the saddle point, 
the value at the saddle point of the free energy and the classical rho mass remain the same. We 
will also require that the residuum of the pion propagator at = is unity, there is no explicit 
term, and the coupling of the X^^ term is g/2. That is we will require the following renormalization 
conditions: 



fs 



-Xl 



dX 



= 0, 

^ref 

dG-\p) 



dp^ 



£ref,P=0 



0, 



9V 



d^^dX 



E, 



9: 



ef 



^ref 



0, 



1, 



Hi 



1, 



where we denoted 



£:,e/ = {T = o, (X,$) = (X,$)^iJ 



(27) 
(28) 



as the reference environment. 



4.2 Leading order 

Here we determine the 0{N^) counterterms, and give the leading order free energy, as well as the 
leading order x propagator. 

The free x propagator is C(iV°), but there are loop corrections of the same order. Since these 
loop corrections are not suppressed by the present series expansion parameter, 1/N , they may appear 
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in any diagram in arbitrary number. That means that in the naive form at each level of perturbation 
theory we have infinitely many diagrams. This phenomenon is clearly must be avoided in a well 
defined perturbation theory, there the radiative corrections must be of lower order in the expansion 
parameter than the leading term. We can achieve this goal with help of our counterterm 6H{p): we 
may tune it in a way that it cancels the 0{N^) radiative corrections to the x propagator, or, in other 
words, we require that the complete self-energy is zero for the x-field: 

2 2 

E^^ip, £) = ~^Iip, £) - SHoip) = ^ 5Ho{p) = £), (29) 

where we explicitly signaled that the computation is performed in a generic S = {T,X, $} environ- 
ment. The / integral reads as 

I{p,S)= I G^{p-q)GM , (30) 

it is the bubble diagram at environment £. The integration symbol at finite temperature corresponds 
to 

Po = 27mT, and we apply an UV as well as an IR cutoff to regularize the integrals. We can then 
calculate the result 

I{p,£) = Ioip)+lT{p), (32) 
where the zero temperature part reads 



1 



, eA^ 7-1 
In — + 7 In — — 

7 + 1 



7^ = 1 + (33) 



where = + Ajp, where Ajji is the IR cutoff. The finite temperature part computed with the 
Bose- Einstein distribution function n{uj) is the following: 



oo 



__y..,„(|±M±M!)„M. (34) 



We may wonder that the counterterm we choose is environment dependent. However, as we 
discussed in Section O the physics, represented by the bare Lagrangian, is insensitive to the actual 
choice of the counterterm, as far as the consistency equation (11 01) is fulfilled. Therefore we must 
choose: ^ 

H{p) = 1 + 6Z^^o - SHoip) = 1 + 6Z^,o + jHp, £) (35) 

for the free x-kernel. Here 5Z^fl is part of the bare Lagrangian, and so it cannot be chosen environ- 
ment dependent. But this is actually unnecessary, since the bubble integral is only overall divergent, 
and so the divergences do not depend on the IR quantities. Taking into account the renormalization 
condition (p7|) we shall choose 

-^In — 

327r2 ^ ml 



^Zxfl = -i^,^^-j, (36) 



then H{p) is a finite quantity with H{—m^^ = 1. 

Some remarks in connection with the leading order result: 
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With the procedure above we defined an optimally improved perturbation theory, where all 
net 0{N^) radiative corrections are missing (the diff'erent contributions cancel each other), and 
so, as a consequence, at each level only finite number of diagrams are generated. Therefore at 
higher orders there is no need to repeat this part of the procedure, only ordinary perturbation 
theory should be used. 

One can realize that at asymptotic momenta the kernel H{p) is logarithmic function of the 
momentum. This ensures that the renormalizability conditions of [20j are satisfied. 

Later we will need the "double asymptotic" form of H, when both the external p as well as 
the internal loop momentum (g) go to infinity at the same rate. This form will be relevant 
for the analysis of overall divergences of higher order diagrams. The result is the same as the 
mn = 0, T = case, since these are IR quantities: 

Hasip) = l + —-^ln^ = —^\n—, L = e^mg, (37) 
SzTT^ p^ p^ 

L is the position of the Landau pole. 

The original bare coupling of the 0(N) model is (cf. (jll) and ([6])) A = ?>gl/ {Z^^Z"^). In the 
leading order it reads 

- _ 3^2 ^ ^ 

^ - g^ , eA^ - A , eA^' ^^^^ 
In — - 1 In 



327?^ 9671^ 

where we defined A = ?)g^ as the renormalized quartic coupling. This is the "nonperturbative 
renormalization" formula for the coupling [21] . 

In order to be able to fix the other leading order counterterms, we shall consider the free energy. 
The classical free energy is proportional to A^, but also the first quantum correction are of the same 
order of magnitude. The interactions ((1 — exp(— S''^*'^))) do not contribute here, just the free pion 
gas. In C^^^^^^^ we also have to take into account those terms which contribute at leading order. 
What we have: 

^ = !^*'^ - !Mx^ + |a-*'^ + - ^^A- + i Tr log + ™i) + */„. (39) 

We use that H + 5H = Z^ = 1 + SZ^, and we introduce an environment independent UV regulator 
to write 

I? = ^ - W + ix^' + ^oX - + 1 Tr log ( P^±^] + 5 /o. (40) 



A^ 2 2 2 ^"2 

The value of the go and 5/o counterterms can be determined by the conditions Taking into 

account the fact that our renormalization conditions were fixed at zero temperature, we find: 



A*^ 777,2 

-A^ + In — + gX^i^ In ^ 



mi emi 



cr I [ d^p , p^ + ml 5^x,o^2 



min 



27 (27r)4 p2 ^ 2 

m\ A2 ^ ^'AT^in, < m. 



1 



327r2 



7„2^2 _ 1^ ^ m| 

2 mi 2 emi 



(41) 
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Note that the divergent parts of go and 6fo are environment-independent, as it must be. As it can 
be easily checked, SZ^^ cancels the divergence coming from the second derivative of the free pion 
contribution. Then /o is a completely finite quantity. 

4.3 Next to leading order 

The next-to-leading order renormalization means that we must determine the part of the coun- 
terterms. Due to the efforts we made at the leading order, this order is the same as any other 
conventional perturbation theory: since the problematic XX self energy due to the pion modes is 
always canceled by the 6H counterterm, we will have only a finite number of diagrams to take into 
account. To determine the counterterms we should consider physical quantities and require that the 
divergences cancel in the perturbative expressions of these quantities. 

First we determine the 6Zi, Sml, 6gi and 6Xi counterterms from the pion self-energy to the first 
order: 

= ^l^ GAP - q)G^M + + K + 59iX + + ^ GM- (42) 

The first term on the right hand side at first look is only overall divergent, there are no divergent 
sub-diagrams, since it is only a one-loop diagram. In the overall divergence all momenta go to infinity, 
and one could think that the divergences cannot depend on the temperature, and the dependence 
on the field background is polynomial: exactly of the form what are the second to fifth terms in the 
above expression. But we must not forget, that here the propagator is not the usual quadratic free 
propagator, and so its temperature dependence can be stronger than an exponential suppression (in 
fact, in the present case it is H{p) ~ T'^/p'^). At the same time the presence of the last term is also 
embarrassing. It is clear that it is of order 1/A^: although formally it is a next-to- next-to-leading 
order contribution, but the number of pions (A^ — 1) lifts up this contribution to the next-to-leading 
order. But what is the role of this diagram? 

If we rewrite the x~ X propagator in the original, not-resummed language, the first contribution 
of fj42|) corresponds to the diagram of Fig{T]/a. But this diagram has a sub divergence, which is shown 




a.) b.) c.) 



Figure 1: The pion self energy as expanded in the language of the non-resummed x ~ X propagator 
(a.). Its potentially divergent subdiagram boxed by dotted line (b.) The diagram necessary to cancel 
this subdivergence (c.) 

by Fig.l/b. To cancel this subdivergence, we clearly want a counterterm diagram as appears in 
Fig.l/c, which is exactly the same as the last term in the pion self-energy. The conclusion is that in 
the original language the sum of the first and the last term is overall divergent. This ensures that the 
SZi,Smi,Sgi and 6Xi counterterms are temperature independent, which is needed for consistency. 
Their value can be determined using the first term in (H^ where all momenta in the original. 
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unresummed representation are asymptotic. In the resummed language it means that we should use 
the asymptotic pion propagator, and the "double asymptotic" form of H{p) described in ( 1371) . 

To read off 6Zi we should consider the dependence of the overall divergent part of the first term 
in (jl2]) (the tt-x bubble integral). flH7|) implies that the x-X propagator behaves as G^^{p) ~ 1/ In(p^). 
So we can write for the overall divergent part of the bubble integral 



d Qi — o ~ / -r^ / de — r— r- -, ~ / -1^ + finite x p\ (43) 

[p — q) inq J Ing J 1 + [p^/q'^) + 2p/q cos J mq 

Therefore there is no divergent one-loop contribution proportional to p"^. This implies 

6Zi = 0. (44) 

For the determination of the cutoff dependence of 6m1 and 6gi we include also the pre-factors of 
(!37|) and write the overall divergent part of the bubble integral at zero external momentum 

A r dqq^ 1 



4 P do 1 

«™;+«9.A-H--/-j^j-^ = finite. (45) 



It means 



and 



= iVlnL/A 



1 1 



21nL/A 21n^L/A 



+ ^^^^^A^' ^^^^ 



(^^Mii, = ^Inln^. (47) 

It is interesting to observe that neither 5m? nor 5g is suppressed by powers of the coupling, both are 
at the same order as im? and respectively. 

To have an expression for the 5\i it is simpler to consider another observable, the 4-point function 
of the scalar fields. It can be expressed at zero external momenta as 

= ^ jGAq)G^AQ)GAQ)G^xiQ) - (48) 

At asymptotic momenta GT,{q) ~ 1/g^ and G^^(g) ~ l/ln(g). Therefore the divergence structure is 
determined by 

dqq' f dlnq ^^^^ 



g^(lng)^ J (Ing)^' 
which is also UV finite expression. That means that 

5Xi = 0. (50) 

We should remark that the convergence of the one-loop contribution to the wave function renormal- 
ization and to the four-point function does not mean that these remain finite at any order. So we 
still have to keep 6Z2, 6X2 and all higher counterterms. 
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For the determination of the counterterms qi,5Z^^i and 5fi we have to write up the effective 
potential to 0{N°). We have 



/i 



N 



5m\ 



2 ^^ 2 24 



5Ai 



H — Tr lo£ 
2 ^ 



+ 







+ 



24 



(51) 



This is the complete expression, which contains the cancellation of all possible sub divergences. To 
understand the role of the different terms, the logic is pretty much the same as in case of the pion 
propagator: the trace- log contribution, although formally only overall divergent, if we expand it with 
respect to the original, non-resummed modes, we find that it is, in fact, an infinite sum of higher order 
diagrams, which contain also subdivergences. This is illustrated in Fig J4.3[ The relevant diagrams 

o o 



O „ D 



a. 







Figure 2: The trace-log contribution, expressed with help of the non-resummed propagators (a.). 
The quadratic subdivergence (b.) and its counterterm (c). The quartic subdivergence (d.) and its 
counterterm (e.). 



necessary to cancel the subdivergences are the first and second terms in the second line in (jSTl) . which 
are automatically generated in the present framework. In the complete expression no subdivergences 
survive, only the overall divergent or finite contributions remain. Since in the overall divergent 
part all momenta go to infinity, there is no way to generate a temperature dependent divergence. 
Therefore all the necessary counterterms are also temperature independent. 

The actual determination of the counterterms is analytically very hard, since it will contain 
diagrams up to three loop order (each derivatives of H{p) with respect to the background fields raise 
the diagrammatic loop number). On the other hand if we use numerical techniques, the expression 
of /i contains only one- loop expressions. Since we know that there exists a consistent, temperature 
independent choice for qi,5Z-^ i and we can perform a zero temperature regularized numerical 
calculation, do the derivations numerically, and determine the correct values for the counterterms. 



5 Conclusions 

In this paper the renormalization of the 0(N) model was discussed. Since the radiative corrections 
can be of the same order of the expansion parameter (now 1/A^) as the leading order result, a loop- 
order based perturbation theory is possible only if we include the leading order quantum effects into 
the quadratic Lagrangian. In IPI technique it can be achieved by integrating out the pionic modes 
|25j . in the 2PI case one applies an external propagator field to represent the non-trivial propagation. 
In the present case we used a non- conventional separation of the bare quadratic Lagrangian into free 
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and counterterm parts, where the free propagator contains all the leading order corrections. In 
order to achieve this we need momentum dependent counterterms, too. Since the so-defined free 
propagator satisfies the conditions of [20], this perturbation theory can be treated with the normal 
renormalization techniques. 

Once we have included the leading radiative corrections into the free propagator, all the rest 
is conventional perturbation theory. We have to take care, however, that the pion modes can still 
lift up certain diagrams to higher level, but now it affects only a finite number of diagrams. If we 
analyze the resulting perturbation theory in the language of the original fields, we can see, that the 
unusual terms correspond to resummed sub divergences of the leading order result. In our calculation 
these terms appear automatically, one does not have to perform a detailed BPHZ renormalization to 
obtain them. 

As a result we established the cutoff dependence of the first nontrivial terms of the counterterms, 
and gave a recipe how numerically can the rest be computed. We gave also the expression for the 
renormalized free energy up to next-to-leading order. 
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